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Abstract 

Analytical  models  are  required  for  the  design  and  development  of 
combustors  whose  applications  range  from  gas  turbine  engines  to  industrial 
boilers.  Recent  advances  in  the  solution  of  the  problem  of  predicting 
the  performance  of  combustors  have  been  based  on  finite  difference  and 
stirred  reactor  methods. 

Three-dimensional  finite  difference  methods  were  pioneered  by 
Professor  Spalding  and  co-workers  and  these  techniques  have  been  successfully 
applied  to  the  prediction  of  gas  turbine  combustor  flow  patterns.  Due  to 
the  limitations  of  present  computers  these  methods  cannot  be  extended  to 
completely  include  fuel  spray  dynamics  and  realistic  chemical  kinetics. 

This  difficulty  is  overcome  by  using  the  computed  flow  patterns  to  define 
a network  of  interconnected  stirred  and  plug-flow  reactors.  The  detailed 
kinetic  scheme  presently  consists  of  13  species  undergoing  18  reactions  to 


I 
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represent  the  combustion  of  hydrocarbon  fuels  such  as  kerosine. 


A model 


of  fuel  evaporation  is  incorporated  fdiich  assumes  the  fuel  spray  to  be 


composed  of  21  x 12  ym  size  ranges,  evaporation  is  calculated  using  a forced 
convection  model.  The  mixing  o£  fuel  vapour  and  air  is  modelled  using  a 
micro-mixing  parameter,  based  on  turbulence  dissipation  rates.  The 

overall  method  thus  combines  3-D  fluid  dynamics,  turbulent  mixing,  evaporation 
and  chemical  kinetics. 


The  model  has  been  verified  by  experiments  which  show  that  flow  velocity 
profiles,  chemical  species  (CO,  NO  etc.)  can  be  successfully  computed. 


Introduction 


Combustion  is  one  of  the  most  difficult  processes  to  model  mathematically 

sinc6  it  generally  involves  the  simultaneous  processes  of  three-dimensional 

two-phase  fluid  dynamics,  turbulent  mixing,  fuel  evaporation,  radiative  and 

convective  heat  transfer,  and  chemical  kinetics.  In  order  to  design 

combustors  based  on  fundamental  principles,  a comprehensive  model  incorporating 
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all  these  factors  is  required  CH*  Unfortunately,  although,  we  can  write 
down  a set  of  governing  differential  equations,  the  solution  of  these 
equations  is  much  beyond  the  capability  of  present  computers  and  a cos^romise 
has  therefore  been  sought.  In  this  study,  the  approach  takes  advantage  of 
the  fact  that  there  is  relatively  weak  coupling  between  the  fluid  dynamic 
aspects  of  the  flow  and  the  chemical  kinetics.  That  is,  the  flow  pattern 
is  largely  defined  by  the  geometry  of  the  combustor  and  is  insensitive  to 
the  detailed  chemical  species  present  in  the  flow  provided  the  approximate 
tenq>erature  rise  is  represented.  The  calculation  is  therefore  divided  into 
two  sections,  each  of  which  is  within  the  capabilities  of  present  coiiq>uters. 

It  is  fair  to  ask  what  are  the  advantages  of  such  a procedure  in  view 
of  the  fact  that  combustors  have  been  designed  and  operated  for  many  years 
without  the  benefit  of  these  methods.  The  answer  is  that  cooAustor  develop- 
ment by  cut-and-try  methods  is  expensive  and  does  not  necessarily  produce  the 
optimum  answer.  Present  pressures  to  minimize  both  fuel  consumption  and 
combustion  generated  pollution  has  lead  to  the  need  to  turn  increasingly  to 
fundamental  considerations,  for  example,  to  help  reduce  nitric  oxide  production. 
A key  factor  in  such  solutions  is  the  effective  fuel  distribution,  and  this 
is  determined  by  the  fuel  preparation  technique.  It  is  therefore  essential 
that  the  model  incorporate  an  accurate  representation  of  the  fuel  droplet 
(or  particle)  size  distribution  and  solves  the  relevant  equations  for 
trajectory  and  evaporation. 

The  first  stage  of  the  calculation  therefore  consists  of  predicting  the 
three-dimensional  flow  pattern,  and  here  we  have  extended  and  applied  the 
finite  difference  methods  pioneered  by  Professor  Spalding  and  co-vorkers  at 
Imperial  College.  The  next  section  of  the  paper  discusses  how  these 
methods  have  been  implemented  and  verified  using  a typical  gas  turbine 
combustor  geometry,  see  Fig.l.  In  order  to  include  the  effect  of  temperatur 
rise  on  the  flow  pattern,  the  model  includes  a crude  two  step  chemical 
kinetic  scheme  and  a crude  droplet  evaporation  scheme.  The  relevant  output 
consists  of  the  distribution  throughout  the  combustor  of  velocity  coiiq>onents , 
turbulent  kinetic  energy,  turbulence  dissipation  rates,  temperature,  pressure 
and  heat  transfer  rates. 


The  second  stage  of  the  calculation  derives  from  chemical  plant 
flowsheet  techniques.  It  consists  of  a procedure  which  represents  the 
combustor  as  a network  of  interconnected  stirred  and  plug  flow  reactors, 
and  includes  a detailed  kinetic  scheme  for  the  chemical  species  to  be 
considered.  A more  detailed  model  of  the  fuel  evaporation  and  mixing 
rate  is  also  built  into  this  part  of  the  computation  since  only  the  fuel 
which  has  evaporated  and  mixed  with  the  air  can  take  part  in  the  chemical 
reaction.  The  objective  of  this  stage  of  the  computation  is  to  predict 
the  coiid>ustion  efficiency  and  pollution  levels  produced  by  the  particular 
combustor  design.  An  important  aspect  of  such  a procedure  is  that  it  should 
be  capable  of  predicting  the  trend  of  dependence.  For  example,  it  is 
important  to  be  able  to  predict  the  effect  which  a finer  fuel  spray  may  be 
expected  to  have  on  the  production  of  nitric  oxide  at  some  particular 
combustor  throughput. 

The  link  between  the  first  and  second  stages  of  the  calculation 
consists  of  identifying  the  appropriate  reactor  network  from  the  computed  flow 
pattern.  The  basic  hypothesis  used  for  this  link  is  that  mixing  at  the 
molecular  level  is  carried  out  by  the  movement  of  molecules  between  turbulent 
eddies  having  different  concentration.  Since  the  molecules  also  carry 
momentum  and  hence  dissipate  the  turbulence,  it  is  assumed  that  mixing 
rate  is  proportional  to  turbulence  dissipation  rate.  The  output  of  the  first 
stage  (flow  pattern)  calculation  gives  the  distribution  of  turbulence  dissipation 
with  high  values  in  the  vicinity  of  the  jets  where  air  is  introduced  into  the 
combustor.  These  regions  of  high  turbulence  dissipation  thus  become  well 
stirred  reactors  (WSR)  of  corresponding  volume,  and  the  degree  of  stirring 
or  mixing  rate  factor  is  determined  from  the  total  dissipation  Cor 
energy  transfer)  within  the  reactor.  Reactors  having  low  mixing  rates  are 
represented  as  plug  flow  reactors  (PFR)  and  it  is  assumed  that  reactions  taking 
place  there  are  confined  to  the  material  which  is  already  mixed.  The  flow 
rates  and  interconnections  between  the  reactors  forming  the  network  are 
determined  directly  from  the  overall  flow  pattern.  The  link  between  the  two 
stages  of  the  computation  is  therefore  complete  and  could  be  automated  in  the 
computer  if  required,  however  we  have  not  yet  made  this  step.  Indeed  the 
purpose  of  the  present  study  is  to  demonstrate  the  validity  of  the  separate 
stages. 
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In  some  cases,  the  region  in  the  Immediate  vicinity  of  the  fuel 
nozzle  is  a major  source  of  pollution;^  since  it  is  unnecessarily  extravagant 
on  computer  time  to  study  this  region  in  detail  using  the  preceding 
procedures,  a separate  axi-'syiimetrical  two  phase  flow  analysis  with  detailed 
droplet  size  distribution  representation  is  used  (2) , however  space 
precludes  the  inclusion  of  this  complementary  study  herein. 


Finite  Difference  Procedure  for  Flow  Pattern  Prediction. 


The  finite  difference  procedure  employed  in  this  work  (3)  uses  the 
cylindrical  polar  system  of  coordinates  to  predict  the  complex  three- 
dimensional  swirling  and  reacting  flow  inside  a combustor. 

The  numerical  scheme  for  solving  the  governing  non-linear  equations  with 
arbitrary  boundary  conditions  coupled  with  mathematical  models  of  turbulence 
and  combustion  has  the  pressure  and  velocities  as  the  main  flow  variables; 
facilities  exist  to  handle  both  physically-controlled  and  kinetically 
controlled  reactions  and  a six-flux  radiation  model. 


The  study  can  undertake  problems  in  Cartesian  and  cylindrical  polar 
coordinates.  Of  these  two  systems  the  cylindrical-polar  coordinates  can 
be  regarded  as  the  more  general  one,  for  the  Cartesian  coordinates  can  be 
derived  from  it  by  siting  the  flow  domain  far  away  from  the  axis.  To 
simplify  the  discussion  all  following  equations  will  be  presented  in 
cylindrical  polar  coordinates  with  the  understanding  that  the  corresponding 
form  of  the  equations  can  be  readily  derived  by  letting  all  the  r's 
approach  infinity  and  replacing  r30  by  3z. 


(a)  Equations  of  continuity  and  momentum. 
Continuity: 


3p 


H + div  (p  v)  - 0 


(1) 


u - momentum: 
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r - momentum: 


♦ div  (p  vu  - ygradv) 


momentum 


+ div  (p  vw-  pgradw) 


(rpu(|>)  + 


div  (ugraditi) 


In  the  above  equations,  u,  v and  w are  the  component  velocities  in  the 
directions  x,  r and  6 respectively;  p is  the  pressure;  p and  p are 
respectively  the  density  and  viscosity  of  the  fluid  mixture. 


For  turbulent  flows  which  are  frequently  encountered  in  cooDbustion  studies 
it  shall  be  assumed  that  the  same  equations,  are  also  valid  provided  that  one 
uses  time-mean  values  for  all  the  flow  variables  and  fluid  properties  and 
that  y is  now  the  effective  viscosity  which  is  the  molecular  viscosity  aug- 
mented by  the  turbulent  contribution.  The  latter  is  derived  at  each  point 
from: 


where  C is  a universal  constant 


The  hydrodynamic  turbulence  model  applied  here  la  a two  equation 
model  of  turbulence  known  as  the  k've  model.  Reference  (4).  It  entails 
the  solution  of  two  transport  equations  for  turbulence  charactiiristics 
namely  that  of  k,  the  local  kinetic  energy  of  the  fluctuating  motion  and 
e,  the  energy  dissipation  rate.  Knowledge  of  k and  e allowa  the  length 
scale  to  be  determined  and  also  the  effective  viscosity  (as  above)  from  which 
the  turbulent  shear  stresses  can  be  calculated. 

The  differential  transport  equations  for  k and  e are  as  follows: 

+ div  (p  vk  - grad  k)  » Gj^  - pe  (6) 

+ div  (p  ve  - r ^effgrad  e)  » ^ 

In  the  above  equations  the  generation  term  for  k,  is  given  bv: 
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r ^ J 

and 

are  constants. 

^k,eff 

exchange-coefficients . 

Recommended  values  for  the  constants  appearing  in  the  above  equations 

are: 

'd  ‘=2 

0.09  1.93  1.92  O) 


(b)  Reaction  Models 

Conservation  equation  for  a chemical  species  j ; 

3(pm»)  ^ 

— + div  (p  vm^  - Tjgrad  mj)  ■ Rj 

where  m^  is  the  mass  fraction  of  the  chemical  species  j»R.j  tbe  mass  rate 
of  creation  of  species  j by  chemical  reaction,  per  unit  volume,  and  Tj  is 
the  exchange  coefficient. 


In  a multi  component  ayartem  aimplifications  can  he  introduced  through 
the  use  of  the  concept  of  a simple  chemically  reacting  system.  This  model 
is  based  on  the  assumptions  that:  the  fuel  and  air  react  chemically  in  a 
unique  proportion;  that  the  effective  diffusivities  of  all  chemical  species 
are  equal;  that  the  reaction  is  a single  step  with  no  intermediate  coiiq>ounds. 


These  suppositions  make  it  possible  to  solve  only  one  equation  for  the 
variable  f defined  as  the  mass  fraction  of  fuel  in  any  form  (burnt  or  unhurnt) 
for  diffusion  flames  and  two  equations  with  variables  m^^  and  f respectively 
for  the  kinetically  influenced  flames. 

The  quantity  f is  defined  mathematically  as: 


f = (v  - V )/(v. 

ox  lu 
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(m-  8 = stoichiometric  constant, 

fu  inlet’ 


(10) 


The  quantities  (m,  ^ and  (m  ).  , ^ are  the  mass  fraction  of  fuel  and 

fu  inlet  ox  inlet 

oxygen  in  the  fuel  and  oxidant  stream  respectively. 

The  mass  fractions  of  fuel  (for  the  diffusion  controlled  case)  and 
oxygen  are  related  to  the  mixture  fraction  according  to: 


0 < f $ f 


St 


m 


fu 


m 


ox 


I > f 5 f 


st 


m 


fu 


(f  - 


m 


ox 


In  the  above  equations  f^^  is  the  stoichiometric  value  of  the  mixture 


fraction  and  is  given  by: 


St 


1/ 


^ ^ inlet 

ox 


(12) 


The  mass  fraction  of  nitrogen  and  the  con^ustion  products  are  treated 
similarly. 

The  influence  of  turbulence  on  reaction  rates  is  taken  into  account  by 
employing  the  eddy  break  up  model  of  Spalding  (5) . 


The  reactloa  rate  in  this  case  is  taken  to  he  the  smaller  of  the  two 
expressions  given  hy  the  familiar  Arrhenius  formulationand  the  eddy  break, 
up  model.  The  latter  is  conveniently  described  as 


where  C is  a constant  and  g represents  the  local  mean  square  concentration 
R 

fluctuations. 


Radiation  Effects 


The  effect  of  radiation  in  the  mathematical  model  are  accounted  for  by 
reference  to  the  six-flux  model  of  radiation.  The  differential  equations 
describing  the  variations  of  the  fluxes  are:  Reference  (6). 


CI+J+K+L+M+N)  ] 


CI+J+K+L+M+N)  ] 


(I+J+K+L+M+N) 


(I+J+K+L+M+N) 
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(I+J+K+L+M+N) 


radiation  flux  in  the  direction  of  positive  r 

radiation  flux  in  the  direction  of  negative  r 

radiation  flux  in  the  direction  of  positive  x 

radiation  flux  in  the  direction  of  negative  x 

radiation  flux  in  the  direction  of  positive  6 

radiation  flux  in  the  direction  of  negative  6 

absorption  coefficient, 
scattering  coefficient. 
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E = aT^  hlaclc  body  emissive  power  at  the  fluid  temperature; 
a S Stefan-Boltzmaim  constant. 

The  composite  fluxes  defined  as: 

= i (I  + J) 

= j (K  + L)  as} 

a j CM  + N) 

are  employed  to  eliminate  I,  J,  K,  L,  M and  N from  the  previous 
equations  to  yield  three  second-order  ordinary  differential  equations. 

(d)  Solution  procedure. 

The  previous  governing  partial  differential  for  mass,  momentum,  energy 
and  species  concentrations  are  elliptic  in  nature  and  can  be  conveniently 
presented  in  the  general  form: 

+ div  (3  t(|)-  rgrad<j)  ) = (16) 


In  the  above  equation  <fi  identifies  the  dependent  variable;  3 is 
identically  equal  to  either  the  mixture  density  p or  zero;  F is  the  appropriate 
exchange  coefficient  for  the  variable  (j>;  and  is  the  source  term. 

Table  (1)  summarizes  the  expressions  for  the  symbols  3,  F and  S^. 

The  equations  are  first  reduced  to  finite  difference  equations  by 
integrating  over  finite  control  volumes (Reference  3)  and  then  solved  by  a 
procedure  described  in  detail  in  Reference  (7)  for  three  dimensional  parabolic 
flows. 

It  will  be  sufficient  for  the  purposes  of  this  paper  to  summarise  the 
main  features  of  the  solution  procedure. 

The  numerical  scheme  is  a semi-implicit  iterative  one  which  starts  from 
given  initial  conditions  for  all  the  variables  and  converges  to  the  correct 
solution  on  the  completion  of  a number  of  iterations. 
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Each  Iteration  performs  the  following  steps: 

i)  The  u,  V and  w momentum  equations  are  solved  sequentially 
with  guessed  pressures. 

ii)  Since  the  velocities  at  this  stage  do  not  satisfy  the  continuity 
equation  locally, a "Poisson-type"  equation  is  derived  from  the 
continuity  equation  and  the  three  linearised  momentum  equations. 
This  "Poisson"  equation  is  then  solved  for  corrections  to  the 
pressure  field  and  the  consequent  corrections  of  the  three 
velocities.  (SIMPLE  algorithm). 

iii)  The  k and  e equations  are  then  solved  using  the  most  up  to  date 
values  of  velocities. 

iv)  The  iteration  is  coiiq>leted  upon  solution  of  the  concentration 
and  all  the  remaining  equations. 

One  point  to  note  here  is  the  implementation  of  the  cyclic  boundary 
condition  in  the  0 direction  compatible  with  the  nature  of  the  swirling  flows 

(e)  Results  of  the  Flow  Analysis  (Computed  and  Experiment) . 

A simple  gas  turbine  combustor.  Fig. 2.,  was  used  to  study  the 
comparison  between  computer  flow  predictions  and  the  experimentally  observed 
values. 

The  air  stream  through  the  swirler  enters  the  chamber  with  finite  swirl 
velocity  while  the  primary,  secondary  and  the  dilution  flows  are  each 
injected  through  6 equally  spaced  circular  orifices  on  the  periphery  giving 
rise  to  the  three-dimensional  nature  of  the  problem. 

There  were  no  wall  cooling  slots  and  the  combustor  was  operated  in  a 
large  plenum  so  that  no  axial  velocity  components  were  included  in  the  radial 
jet  flows.  The  finite  difference  grid  network  reported  here  consisted  of  27 
transverse  planes,  ^ radial  planes  and  18  circumferential  surfaces  (Fig. 2.). 

The  cyclic  boundary  condition  mentioned  above,  can  be  employed  to 
consider  only  a 60°  sector  with  the  particular  input  data  given  below: 


Total  air  flow  rate  - 0.1275  kg/aec. 

Swirler; 

Air  flow  rate  « 0.0016575  kg/sec  for  the  60**  sector. 

Fuel  flow  rate  - 0.0  (cold),  0.000178  (hot) 

Swirl  - 0.8 

Inlet  temp.  ■ 351  k 

Axial  velocity  = 27.6  m/s 

Swirl  velocity  » 14  22  m/s  (Forced  vortex,  w 'Vi  radius) 

Injection  Details; 

Primary  Injection. 

Flow  rate  = 0.005421  kg/sec.  per  orifice. 

Injection  velocity  = 138.90  m/sec. 

Secondary  Injection. 

Flow  rate  « 0.006355  kg/sec.  per  orifice. 

Injection  velocity  »=  139.1  m/sec. 

_0ilv^^iu 

Flow  rate  «■  0.007820  kg/sec.  per  orifice. 

Injection  velocity  =»  138.7  m/sec. 

Injection  temperature  (for  all  ports)  » 351  K. 

It  is  difficult  to  represent  Che  full  detail  of  the  flow  in  a paper 
such  as  this,  therefore  a few  profiles  have  been  selected  to  indicate  the 
character  of  the  results. 

In  Fig. 3,  the  calculated  radial  velocities  are  shown  at  two  radii 
for  a longitudinal  cross-section.  The  decay  of  the  radial  jet  velocity  and 
the  small  radial  outflow  between  the  holes  will  be  noted. 

In  Fig. 4.  the  axial  velocities  are  shown  for  a number  of  axial  locations. 
As  in  the  case  of  Fig.  3,  this  graph  is  for  a section  through  the  centre  of 
the  holes,  and  the  radial  expansion  of  the  jets  is  clear.  Particularly 
noteworthy  is  the  negative  fl.  oundary  in  the  primary  zone,  which  indicates 
the  boundary  of  the  recirculation  zone.  This  location  of  the  recirculation 
zone  was  simply  verified  on  the  test  rig  by  introducing  a sodium  tracer  into 
the  flame.  When  the  tracer  was  introduced  into  the  recirculation  zone,  the 
whole  exit  flow  was  coloured  yellow,  whereas  only  a small  region  was  yellow 
when  it  was  introduced  outside  the  recirculation  zone.  The  axial  velocity 
profile  at  the  exit  will  be  discussed  in  more  detail  below. 
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Fig.S.  shows  the  radial  velocity  profile  at  three  radii  plotted  at  a cross- 
section  through  the  primary  jeta.  The  slight  skewness  caused  by  the  swirl  can 
be  detected,  and  the  swirl  velocity  profile  ia  also  plotted  on  the  same  figure. 

In  Fig. 6 the  turbulence  dissipation  contours  are  plotted  for  an  axial  cross- 
section.  Two  sets  of  contours  are  shown  corresponding  to  sections  through  the 
jets  and  at  the  edge  of  the  jets  respectively.  This  figure  clearly  illustrates 
the  fact  that  dissipation  Cand  mixing)  is  concentrated  in  the  vicinity  of  the  jets, 
and  corresponds  to  the  stirred  reactor  locations. 

It  is  difficult  to  make  accurate  three  dimensional  velocity  measurements  in  a 
small  combustor.  The  problem  is  that  a 5 or  7 hole  spherical  probe  occupies 
too  much  space  in  the  strong  shear  flow  within  the  combustor.  It  was  therefore 
decided  that  the  best  comparison  location  between  the  predicted  and  measured  flow 
fields  would  be  obtained  at  the  combustor  exit.  In  Fig. 7a,  these  experimental 
and  theoretical  exit  axial  flow  profiles  are  plotted  and  it  can  be  seen  that  the 
agreement  is  remarkably  good.  In  Fig. 7b,  the  more  con^lex  factor  of  exit 
turbulence  intensity  is  plotted  and  again  experimental  and  theoretical  values 
are  compared.  The  experimental  values  were  obtained  with  a photon-correlation 
laser  doppler  anemometer.  Again  the  agreement  is  remarkably  good  considering 
the  difficulty  of  the  experiment  and  the  complexity  of  the  prediction,  however 
the  results  suggest  that  the  model  slightly  overestimates  the  turbulence.  In 
Fig. 7c  the  exit  temperature  profiles  are  compared  and  again  the  model  successfully 
predicts  the  poor  exit  traverse  with  a temperature  minimum  close  to  the  axis. 

Stirred  Reactor  Modelling. 

As  pointed  out  previously,  the  second  stage  of  the  overall  cosohustor 
modelling  procedure  is  to  riodel  the  flowfield  by  an  appropriate  network  of  stirred 
reactors.  Consequently  a sub-model  is  required  for  a Well  Stirred  Reactor,  WSR, 

(N.B.  not  perfectly  stirred)  which  includes  the  internal  processes  of  fuel  evap- 
oration, mixing  and  complex  chemical  kinetics.  It  is  possible  to  use  "global" 
reaction  rate  data  in  such  a system  to  predict  combustion  performance  and  combustor 
stability  (ref. 8),  however,  if  prediction  of  pollutant  emissions  is  to  be  attempted 
then  detailed  chemical  kinetic  schemes  are  required.  In  addition,  to  model  tbone 
sections  of  the  flowfield  where  no  mixing  is  taking  place,  a Plug  Flow  Reactor 
model  is  required,  this  can  be  achieved  by  a sequence  of  differentially  small 
W.S.R.s. 

Previously  stirred  reactor  models  of  combustors  have  been  mainly  confined 
to  homogeneous  combustion  with  "global"  reaction  kinetics  (ref. 9),  although 
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Gaseoua  phase  mass  balance  - 02  * (Wdt 

Liquid  phase  mass  balance  - ^2  **  FE  ■ de/dt 

Mass  balance  on  the  unmixed  fluid  in  the  gaseoua  phase 
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T^/tp  (unmixedness  parameter) 


Thus  for  a given  FE,  tgp  and  feed  conditions  (.25)  defines  the  proportion  of 
the  steady  state  reactor  gas  phase  which  is  unmixed.  Therefore  the  reactor 
composition  prior  to  reaction  may  be  expressed  as  follows:- 


encouraging  results  were  obtained  within  the  limitations  of  this 
approach  (ref. 8).  Evaporation  effects  were  later  included  with  some 
success  (ref. 10). 
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Although  the  three  processes  of  evaporation,  mixing  and  reaction  must 
occur  simultaneously  in  a combustor  it  is  convenient  to  consider  that  in 
a steady  state  condition  these  processes  occur  in  series.  Thus  the  fuel 
must  first  evaporate,  then  mix  and  finally  react.  This  approach  allows  us 
to  calculate  the  reactor  gaseous  phase  composition  after  evaporation  and 
mixing  have  taken  place,  this  constitutes  the  homogeneous  feed  to  the 
reactor.  A general  combustion  scheme  which  summarizes  this  is  shown  in 
table  1. 

1 

The  feedstream  to,  or  product  stream  from  any  reactor  is  assumed  to  be 
composed  of  any  or  all  of  the  eight  species  in  the  above  scheme.  Figure  8 
shows  the  composition  of  the  general  two  phase  steady  state  reactor  with 
the  liquid  phase  shown  coalesced  for  convenience.  Transfer  from  the  liquid 
phase  to  the  gas  phase  is  represented  by  the  mean  fuel  evaporation  rate, 

Fk.  In  addition  it  is  assumed  that  the  mean  residence  time  of  each  phase 
in  the  reactor  is  the  same.  This  assumption  is  incompatible  with  the  exist- 
ence of  a relative  velocity  between  the  gas  and  the  fuel  droplets  but 
greatly  simplifies  the  analysis  and  calculation  of  fuel  distribution  around 
any  particular  reactor  network.  It  is  not  essential  to  assume  this  however, 
and  the  analysis  could  be  modified  to  incorporate  unequal  phase  residence 
times.  Transfer  of  fluid  from  the  unmixed  state  to  the  isixed  state  is 
assumed  to  take  place  at  the  following  rate:- 

(mass  of  unmixed  fluid) /t^ 

where  is  the  characteristic  turbulence  dissipation  time. 

Total  reactor  mass  =m+e  (17) 


The  concentrations  of  the  untnixed  species  can  be  obtained  by  performing 
the  relevant  species  mass  balances.  For  example,  unmixed  fuel  vapour  mass 
balance  is  as  follows 
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The  expressions  for  the  intermediate  concentrations  of  the  mixed 
reactants  and  combustion  products  are  obtained  by  performing  similar  mass 
balances  on  the  mixed  portion  of  the  reactor  gas  phase.  For  example,  the 
mixed  fuel  vapour  mass  balance  is  as  follows 
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Using  (25)  and  (28) 
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Having  defined  the  intermediate  composition  of  the  gaseous  phase,  i.e. 
after  mixing  and  evaporation,  tvo  further  balances  must  be  made  to  determine 
the  final  reactor  composition.  These  are  the  species  and  energy  balances 

for  the  PSR  for  the  steady  state.  During  the  reaction  stage  the  mixed 

* 

gaseous  phase  concentrations,  y , are  transformed  to  the  final  concentration 
Y;  the  unmixed  gas  phase  concentrations,  u,  do  not  change  of  course  although 
they  do  contribute  to  the  physical  properties  of  enthalpy,  specific  heat  and 
density. 

Chemical  reaction,  mixed  species  mass  balance:- 


KTvTf:  - ♦„>  * ^ 


- 0 


(31) 


where  i « 1,  Ml,  > total  nunfter  of  mixed  gaseous  species. 
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Chemical  reaction,  gaa  phase  energy  balance 


0 for  adiabatic  operation 


The  species  kinetic  production  term  is  given  by 


The  forward  and  backward  reaction  rates  are  related  to  the  reactor  gas  phase 
species  concentrations  by  the  following:- 


where  X.  is  a third  body  in  a dissociation  reaction 


The  forward  reaction  rate  constants  are 


The  backward  reaction  rate  constants  are  fixed  by  the  equilibrium  constants: 


The  system  of  equations  is  conq>leted  with  the  equation  of  state :■> 
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The  final  reactor  gas  phase  composition  has  thus  been  defined  and  the  final 
overall  concentrations  may  be  expressed  as  follows 

‘^f  ^ '1'^  Uj  + (1  - <l>y)  Yf  (41) 

and  similarly  for  the  other  species. 

The  complete  set  of  equations  characterises  a heterogeneous  stirred  reactor 
in  steady  state  operation,  they  reduce  to  the  equations  for  the  homogeneous 
case  if  B is  set  to  zero.  These  equations  are  very  non-iinear  due  to  the 
exponential  dependence  of  the  reaction  rates  in  temperature,  additionally  FE 
is  a complex  function  of  temperature  and  staytime,  t^,  hence  the  WSR 
equations  must  be  solved  iteratively.  The  technique  used  is  based  on  the 
numerical  method  PSR  solution  developed  by  Osgerby  (ref. 11),  in  which  a 
Newton-Raphson  correction  procedure  is  employed  to  converge  onto  the  PSR 
solution,  from  an  initial  guess.  The  data  required  for  the  solution  are  a 
kinetic  scheme,  rate  data,  thermodynamic  data,  feed  conditions  and  an 
initial  guess.  The  initial  guess  is  provided  by  an  equilibrium  calculation 
and  the  kinetic  scheme  used  is  shown  in  table  2. 

In  order  to  calculate  fuel  evaporation  rates  a timestep  technique  was 
devised,  the  spray  is  assumed  to  consist  of  21  size  intervals  each  of  12  y. 


I , each  of  which  is  represented  by  the  interval  mean  diameter.  The  percentage 

I of  the  fuel  spray  remaining  unevaporated  at  any  elapsed  time  decreases  non- 

'•  linearly  with  time,  therefore  the  evaporation  rate  of  the  spray  is  not 

I I constant  through  the  reactor.  The  spray  mean  evaporation  rate  is  obtained 

i by  calculating  the  total  fuel  evaporated  during  the  staytime  in  the  reactor 

! and  dividing  by  the  staytime.  Initially  the  droplets  have  a velocity  relative 


to  the  gas  straaa,  howaver  as  tha  drag  forces  acting  on  the  droplets  are 
inversely  proportional  to  the  diameter  the  small  droplets  rapidly  assume  the 
local  gas  velocity  whereas  the  larger  droplets  tend  to  retain  their  own 
velocity.  There  are  two  modes  of  combustion  possible  for  an  evaporating 
fuel  droplet;  droplet  diffusion  flames  and  droplet  wake  flames.  The  wake 
flames  are  generally  blue  due  to  good  mixing  prior  to  combustion  whereas 
diffusion  flames  are  typically  yellow  due  to  soot  formation.  The  velocity 
necessary  to  cause  a transition  from  diffusional  to  wake  burning  is  a 
strong  function  of  the  local  oxygen  concentration  and  falls  to  zero  at  oxygen 
concentrations  in  the  range  14  - 16Z,  thus  at  such  oxygen  levels  a diffusion 
flame  cannot  exist  if  there  is  any  relative  velocity.  This  is  generally  the 
case  in  gas-turbine  combustion.  Using  the  evaporation  model  of  Wise  et 
al  (ref. 12)  we  can  derive  the  static  evaporation  rate: 

4 TT  Ar 

Aj.  = _ In  (1  + (42) 

C 

P 

where  B = C (T  - TJ/L 

ev  p ' CO  L 

2/3 

and  A - 1.432  x lo"^  Cp  (T  - 44.67)  W m"^  k“^ 


To  allow  for  the  effects  of  droplet  dynamics  an  enqpirical  correlation  of  the 
type  suggested  by  FrBssling  is  used  (13):- 


where 


- Ag  (1  + 0.244  Re^) 


Re  = 2 v^^j  r Pg 
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In  order  to  incorporate  drag  effects  into  the  model  an  expression  is  required 
for  the  acceleration  experienced  by  an  individual  droplet,  the  expression 
derived  by  Vincent  (14)  is  used:- 
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The  only  information  needed  to  allow,  calculation  of  the  evaporation  rate 
is  the  initial  droplet  size  distribution  and  the  droplet  initial  velocity, 
these  are  normally  provided  by  correlations  derived  experimentally  for  the 
atomizer  in  use.  In  our  studies,  we  use  the  laser  diffraction  dropsize 
distribution  meter  which  we  have  developed,  to  characterize  the  spray 
accurately  CIS) . 

Before  vaporized  fuel  and  oxidant  can  react,  the  respective  molecules 
must  be  brought  into  intimate  contact,  the  physical  processes  involved  are 
termed  mixing.  Mixing  is  important  under  condiustion  conditions  since  it 
is  usually  the  rate  determining  step,  however  it  is  the  most  difficult 
process  to  model  mathematically.  The  principle  source  of  mixing  energy  in 
a gas  turbine  combustor  is  the  pressure  loss  across  the  turbulence  generator, 
the  can.  Since  the  velocity  and  concentration  fluctuations  decay 
simultaneously  it  is  proposed  that  the  degree  of  mixing  is  equal  to  the  degree 
of  turbulence  dissipation  within  the  flow  system.  An  energy  balance  is 
performed;- 


Pressure  drop 
across  baffle  * 

dP/q  - ^av^*^ 


Energy  "held" 
in  flow  + 

velocity  profile 


turbulence 
kinetic  energy 


dissipation 
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It  is  proposed  that  a characteristic  dissipation  time  t^/u*^^, 

where  C = constant  (unity),  ■ mean  size  of  energy  containing  eddies 

( 0.2  X),  u*  • maximum  value  of  the  r.m.s.  velocity  fluctuations. 

' _ max 

But  T - X/U  (X  “ 10  X)  .’.  T__  = 50  (u'/O)  , thus  using  the  energy 

s bu  I max 

balance  (45),  this  yields  •=  50  (AP/3q)  . 


Thus  the  unmixedness  parameter,  Xg^,  used  from  equation  25  onwards  can 
be  related  to  the  system  geometry.  This  parameter  has  high  values  for  good 
mixing  (>  200)  and  the  values  (<  5)  cause  "blowout"  due  to  inadequate  mixing 
(fig. 9).  A mathematical  model  of  each  process  has  now  been  described  and 
these  are  combined  using  the  equations  derived  earlier  (17  41).  The 

equations  are  solved  in  the  following  way,  firstly  a gas  temperature  is 
estimated  so  that  an  initial  evaporation  rate  for  the  reactor  can  be 
calculated.  Then  using  the  above  equations  the  homogeneous  feed  to  a 
reactor  can  be  calculated  and  an  equilibrium  calculation  performed  to 
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generate  the  etarting  values  for  the  iterative  calculation.  Several 
iterations of  the  chemical  species  equations  are  performed  until  the  mass 
convergence  test  is  satisfied,  then  the  energy  convergence  test  is  applied, 
If  this  is  not  satisfied  then  the  tenqperature  is  corrected  by  (H  ■'  H)/C^ 
followed  by  recalculation  of  the  fuel  evaporation  rate  and  homogeneous 
feed  condition,  after  which  the  Newton-Raphson  scheme  for  solving  the 
chemical  species  equations  is  re-entered.  It  was  found  that  if  the  fuel 
evaporation  rate  was  recalculated  every  time  the  teiq>erature  was  corrected 
then  a slight  instability  is  introduced  into  the  iteration,  therefore  the 
evaporation  rate  is  only  calculated  for  the  first  ten  iterations  by  which 
time  the  correction  is  less  than  5K. 

Ue  have  now  derived  a method  for  solving  an  individual  reactor;  we 
have  to  devise  a sequence  of  reactors  to  represent  the  combustor  under 
consideration.  There  are  several  methods  which  can  be  used  to  determine 
these  sequences :-alQualitative  observations  of  combustor  performance,  e.g. 


for  a typical  gas  turbine  combustor  this  leads  to  a series  model 
consisting  essentially  of  a stirred  reactor  followed  by  a plug  flow  reactor. 
b)The  distribution  of  "macro  scale  mixing",  derived  from  theoretical  and 

experimental  tracer  response  fractions  (e.g.  water  modelling  and  Argon  tracer 
studies). c)  The  distribution  of  "micro  mixing"  energy  within  the  flow;  in 
this  method  the  stirred  reactors  are  placed  in  the  flow  regions  where  high 
levels  of  turbulence  energy  exist.  As  already  discussed  above,  the  3-D 
finite  difference  procedure  is  used  to  determine  flow  patterns  and  turbulent 
energy  distribution  and  hence  derive  the  volumes,  flow  rates  and  inter- 
connections of  the  appropriate  stirred  reactor  network.  An  example  of  a 
reactor  model  of  a gas  turbine  combustor  is  shown  in  fig. 10. 

As  a simple  illustration  of  the  stirred  reactor  network  approach  we 
will  discuss  its  application  to  a novel  design  of  "blue  flame",  low  polluLtnn 
burner  built  at  Sheffield  (fig. 11).  Air  is  added  via  a Coanda  ejector, 
which  causes  controlled  recirculation  of  combustion  products  thus  reducing 
the  oxygen  concentration  in  the  vicinity  of  the  spray  leading  to  wake 
burning.  The  isothermal  burner  flow  pattern  was  determined  by  hot  wire 
anemometry  measurements.  (fig. 12).  The  general  flow  pattern  under  burning 
conditions  was  assumed  to  be  the  same  as  that  under  isothermal  conditions 
although  flowrates  were  corrected  for  density  effects.  Since  the  mean 
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residence  time  in  the  Coanda  throat  la  small,  no  stable  condiustion  can 
take  place  there,  and  the  throat  ia  assumed  to  be  an  isothermal  mixer  of 
atomized  spray,  feed  air  and  recirculated  gases.  WSRl  is  the  main 
flame  zone,  its  volume  is  set  equal  to  that  of  the  truncated  cone 

3 

bounded  by  the  dotted  line  (240  cm  ) . The  unmixedness  parameter, 
for  this  reactor  was  estimated  to  be  300,  i.e.  virtually  perfect  mixing. 

The  secondary  flame  zone  is  represented  by  a poorly  mixed  WSR2  (volume 

3 

280  cm  ),  this  WSR2  was  estimated  to  have  a T--  of  10,  (fig.  13). 

A cooling  of  the  recirculation  flow  around  the  narrow  coanda  unit  annulus  is 
to  be  expected  and  was  apparent  from  temperature  measurements  carried  out. 
Consequently  a heat  exchanger  unit  was  incorporated  into  the  recycle  path 
to  take  account  of  this  (Fig. 14).  The  kinetic  scheme  used  for  the 

calculations  was  that  referred  to  earlier,  table  2,  the  rate  of  reaction  1 
is  given  by:- 


d[C^, 

dt 


5.52  X 10®. exp (- 12900 /T).  [C^2  ^4^^*  ^02^*  (50) 


In  the  experimental  system  water  was  condensed  in  the  sampling  system 
therefore  for  comparison  purposes  the  predicted  water  concentration  was 
reduced  by  a factor  of  10  and  all  other  concentrations  adjusted  accordingly. 
The  initial  size  distribution  was  described  by  the  Rosin-Rammler  expression 
(fig. 15),  the  -Rammler  parameters  were  related  to  the  fuel  injection 

pressure  and  atomizer  fuel  flow  number  using  correlations  of  the  type 
suggested  by  Bowen  and  Joyce  (ref. 16)  based  on  our  own  experimental  results. 
The  spray  initial  velocity  was  calculated  using  discharge  coefficients  based 
on  Tipler's  results  (ref. 17),  assuming  that  all  the  droplets  were  projected 
with  the  same  initial  velocity. 

Results  of  the  Stirred  Reactor  Analysis  (Computed  and  Experiment) . 

The  predicted  values  of  NO,  CO  and  temperature  at  the  burner  exit 
(after  WSR2)  are  plotted  against  airflow  rate  (and  hence  ^ ) for  a given 
set  of  operating  conditions  (Figs. 16, 17, 18, table  3)  with  the  corresponding 
experimental  values.  For  each  of  these  variables  the  correct  trends  and 
order  of  magnitude  is  predicted  and  in  the  case  of  NO  and  CO  the  agreement 
is  particularly  good.  Work  is  still  in  progress  to  extend  the  computed 
results  to  rich  mixtures  and  these  will  be  reported  in  due  course.  The 
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temperature  predicted  is  generally  larger  than  the  measured  temperature, 
by  up  to  200  K,  this  is  due  to  the  measured  temperature  (by  a simple  bare 
thermocouple)  not  having  been  corrected  for  heat  losses.  The  effect  of  fuel 
pressure  on  NO  emissions  is  shown  in  fig. 19  (table  4),  again  very  good  agree- 
ment for  both  trend  and  magnitude  was  obtained. 

Discussion. 

The  major  effort  in  these  studies  has  been  to  reduce  the  computer  time 
, to  manageable  proportions  while  retaining  realistic  combustor  configurations, 
and  in  this  endeavour  we  have  been  relatively  successful.  In  the  finite 
difference  procedure,  an  obvious  variable  is  the  number  of  nodes  on  the  grid 
representing  the  flow  field.  Preliminary  studies  showed  that  a grid  much 
coarser  than  that  used  here  was  inadequate  to  represent  the  gas  turbine 
combustor  flow  field,  and  similar  comments  would  apply,  for  example,  to  an 
industrial  multi-orifice  gas-fired  tunnel  burner.  However  the  system  could 
be  simplified  considerably  if  the  flow  were  axi-symmetric,  as  then  only  two 
dimensions  appear  in  the  analysis. 

A major  weakness  of  the  finite  difference  procedure  is  the  mathematical 
representation  of  the  turbulence  as  isotropic,  whereas  in  the  important  jet 
flow  region  it  is  undoubtedly  anisotropic.  The  significance  of  this  fact 
is  difficult  to  ascertain  since  it  is  probably  compensated  by  the  values  of 
the  coefficients  used  in  the  analysis.  The  fact  that  these  same  coefficients 
give  analytical  results  in  reasonable  agreement  with  a wide  range  of 
experiments  gives  confidence  in  their  use  in  the  systems  described  here. 

It  is  anticipated  that  even  more  sophisticated  turbulence  models  can  be  used 
in  the  same  basic  procedure,  as  they  become  available. 

The  finite  difference  procedure  typically  evaluates  about  twenty  partial 
differential  equations  giving  the  values  of  about  twenty  variables  at  each 
of  3402  grid  points.  The  computer  time  (ICL  1906S)  for  convergence  is  about 
2 hours  for  each  case.  Separate  studies  have  shown  that  the  fuel  spray 
must  be  divided  into  more  than  about  20  size  ranges  if  it  is  to  be  accurately 
represented  in  calculations,  and  if  incorporated  into  the  finite 
difference  procedure  to  this  precision  it  would  almost  double  the  number  of 
equations  to  be  solved  at  each  grid  point,  as  well  as  dramatically 
increasing  the  storage  requirement.  The  computer  time  would  then  increase 
to  the  point  where  the  costs  could  not  be  justified  even  if  the  computer 
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power  and  time  were  available.  Similar  considerations  apply  to  the 
inclusion  of  a full  chemical  kinetic  scheme  into  the  finite  difference 
flow  field  analysis,  since  it  has  been  found  that  about  20  reactions  is 
the  minimum  which  will  adequately  represent  the  combustion  process  if 
pollutants  are  to  be  predicted.  We  therefore  conclude  that  the 
overall  procedure  proposed  here  is  about  the  simplest  that  can  adequately 
represent  all  the  important  physical  and  chemical  factors  controlling 
the  behaviour  of  a real  combustor. 

The  most  important  factor  yet  to  be  determined  in  the  stirred  reactor 
model  is  the  coefficient  of  the  mixedness  factor  At  present,  this 

coefficient  is  assumed  to  be  unity,  however  there  is  hope  that  it  can  be 
determined  more  accurately  by  experiments  in  which  the  probability  density 
distribution  of  the  concentration  fluctuations  are  measured.  An 
appropriate  simplified  analysis  is  given  by  Pope  in  Reference  18. 

Fortunately  the  performance  of  most  combustors  is  not  very  sensitive  to  this 
coefficient. 

A further  weakness  of  the  stirred  reactor  program  is  that  it  cannot,  at 
present,  predict  the  formation  and  agglomeration  of  soot.  This  problem  can 
be  approached  from  two  directions.  In  the  first,  we  assume  that  soot  is 
an  undesirable  product  from  gas  turbine  combustors,  for  instance,  and  the 
foregoing  design  procedure  may  therefore  be  used  to  avoid  the  local  rich 
conditions  under  which  soot  is  formed.  In  the  second,  more  research  is 
needed  to  determine  the  appropriate  rate  equations  for  soot  formation,  and 
these  equations  may  then  be  incorporated  into  the  existing  computer  program 
with  a simple  change  in  the  input  data.  If  the  fuel  to  be  employed  were 
residual  oil  or  coal  which  burn  with  a more  complex  process  than  simple 
evaporation,  then  the  appropriate  module  of  the  computer  program  would  have 
to  be  modified.  However  the  modular  nature  of  the  program  and  simple 
physical  form  of  the  variables  allow  such  changes  to  be  made  relatively 
easily. 


Before  all  the  techniques  reported  above  can  be  applied  with  confidence, 
much  more  comparison  with  experiment  is  required,  perhaps  leading  to  a 
refinement  in  the  procedure.  This  study  thus  reports  a phase  of  progress 
on  a steadily  evolving  subject. 


Conclusions 


1.  The  problem  of  mathematically  modelling  combustion  systems  has  been 
approached  by  a combination  of  finite  difference  and  stirred  reactor  net“ 
work  methods. 

2.  The  finite  difference  model  predicts  the  flow  and  turbulence  field 
and  provides  the  data  for  setting  up  the  stirred  reactor  network. 

3.  The  stirred  reactor  network  model  incorporates  unmixedness  and  fuel 
evaporation  and  predicts  the  combustion  efficiency  and  pollutant  production. 

4.  Experiments  so  far  carried  out  have  verified  that  the  flow  profiles 
and  chemical  sp-tcies  (CO,  NO  etc.)  can  be  successfully  predicted. 

5.  Further  work  is  needed  to  verify  the  method  over  a wider  range  of 
combustors  and  operating  conditions. 

6.  This  study  presents  a phase  of  progress  in  a complex  but  steadily 
evolving  subject. 
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NOMENCLATURE 


i)  Finite  difference  equations 


S’  S’  S 


E 

g 

s 


I,J,K,L,M,N 

k 

I 


m. 


• m » m 
fu’  ox  pr 


q’'.  q’'.  q* 


R, 


rS  r^.  r" 


r 

s 

S, 


T 

t 

u,  V,  w 


Absorption  coefficient  for  radiation. 

Constants  in  the  turbulence  model. 

Constant  in  the  eddy  break-up  model. 

Black  body  emissive  power. 

Mean  square  concentration  fluctuation  of  fuel  species. 

quantity  in  the  generation  term  for  k. 

radiation  fluxes  in  the  r,  x and  6 directions. 

kinetic  energy  of  turbulence. 

length  scale  of  turbulence. 

mass  fraction  of  a species  J. 

mass  fractions  of  fuel,  oxygen  and  product  respectively 
pressure. 

net  radiative  transfer  in  the  coordinate  directions, 
mass  rate  of  creation  of  species  J by  chemical  reaction 


a dependent  variable  for  radiation  fluxes  in  the  three 
coordinate  directions. 


X.  z 


distance  from  axis  of  symmetry, 
stoichiometric  mass  ratio, 
scattering  coefficient  for  radiation, 
source  term, 
absolute  temperature, 
time. 

velocity  components  in  the  x,  r and  z (or  6)  directions 
velocity  vector, 
coordinate  distances. 


mean  kinetic  energy. 

total  number  of  muced  gaseous  species. 

mass  of  reactor  gas  phase. 

total  number  of  gaseous  species. 

temperature  exponent  in  modified  Arrhenius  equation, 
pressure . 

rate  of  production  of  species  i by  chemical  reaction, 
dynamic  head.^ 
uiiiveral  gas  constant. 

Reynold's  number  of  droplet, 
droplet  radius, 
temperature . 

temperature  of  droplet  surface, 
mean  velocity, 
fluctuating  velocity, 
reactor  volume, 
velocity  relative  to  droplet, 
molecular  weight  of  species  i. 

third  body  concentration  in  a dissociation  reaction  j . 

stoichiometric  coefficient  of  species  i as  a product  in 
reaction  j . 

FE/m^  non-dimensionalized  evaporation  rate. 

mass  fraction  of  species  i in  reactor  mixed  gas  phase. 

stoichiometric  coefficient  of  species  i as  a reactant  in 
reaction  j . 

mass  of  reactor  liquid  phase. 

mean  thermal  conductivity  of  gas  phase. 

viscosity. 
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r 

diffusion  coefficient. 

^J.eff 

c 

effective  exchange  coefficient  for  J. 

dissipation  rate  of  turbulence. 

y 

viscosity. 

p 

density. 

0 

Stefan-Boltzman  constant. 

* 

general  dependent  variable. 

e 

coordinate  in  the  cylindrical-polar  system 

ii)  Stirred 

Reactor  equations. 

Pre-exponential  factor  of  the  jth  reaction. 

Reverse  reaction  rate  of  the  jth  reaction. 

rate  constant  of  the  reverse  of  reaction  j . 

transfer  number  in  evaporation  rate  equation  [C  (T-T  )/L] 

P 

overall  gas  phase  mass  fraction  of  ith  species, 
mean  heat  capacity  of  gas  phase, 
droplet  drag  coefficient, 
dissipation  energy  term. 

third  body  efficiency  of  species  i in  reaction  j (usually 
assumed  to  he  1) . 

activation  energy  of  the  jth  reaction, 
fuel  evaporation  rate, 
forward  reaction  rate  of  reaction  j . 
standard  molar  free  energy  of  species  i. 
forward  rate  constant  of  reaction  j . 
rate  of  enthalpy  loss  from  the  reactor, 
specific  enthalpy  of  species  i. 
equilibrium  constant  of  jth  reaction. 


I 

I 


1 

i 

I 


i 


i 

I 


1 

I 

! 

i 


u 


u. 

1 


Superscripts 


index  to  indicate  whether  the  third  body  participates, 
characteristic  dissipation  time, 
reactor  stay  tite. 

unmixedness  parameter, 
proportion  of  gas  phase  which  ia  unmixed, 
mass  fraction  of  species  i in  reactor  unmixed  gas  phase 


1 

* 

Subscripts 

1 

2 

G 

L 

f 

^2 

cs 

i 

j 

E 

E,F 


feed  conditions. 

intermediate  value  Ci.e.  after  mixing  and  evaporation) . 


reactor  inlet, 
reactor  outlet, 
gas  phase, 
liquid  phase, 
fuel . 
nitrogen, 
oxygen. 

general  combustion  species, 
chemical  species  identifier, 
chemical  reaction  identifier, 
evaporation, 
forced  evaporation. 
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GRID  ARRANGEMENT  27  x IB  x 7 x 3 402  TOTAL  POINTS 


Figure  2.  Finite  difference  i:rld  lor  the  combustor. 


RADIAL  PROFILE 

K s A PLANE  THROUGH  HOLES 


Figure  3.  I.ycoming  combustor,  cold  flow  radial  velocity  nrofl.es 
^^:ial  cross-section. 


Figure  4.  Lycoming  combustor,  cold  flow  radir.l  velocity  profiles 
axial  cross-sectior . 
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Figure  5.  Lycoming  combustor,  cold  flow  radial  and  1 argent ial 
velocity  profiles,  transverse  cross-section. 
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Figure  7c 


Figure  7a  Lycoming  combustor,  colfl  exit  velocity  profiles  - predicted  an.' 
measured . 

Figure  7b  Lycoming  combustor,  exit  turbulence  Intensity  profiles  with  combust lo,., 
predicted  and  measured  with  photon  correlation  anemometer. 

Figure  7c  Lycoming  combustor,  exit  temoerature  profiles  with  combustion,  nrertic  e. 
and  measured. 
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Figure  14.  Blue  flame  burner  stlrrec’  reactor  network. 


Figure  15.  ’Measured  atomizer  initial  rtrop.'jize  distribution  charnc  i e^'ist ' cs 
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Figure  16.  Burner  exit  Nitric  Oxide  concentrations,  Measure!  r.nd 
predicted,  vs  air  flow  rate. 
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Recent  advances  in  the  solution  of  the  problem  of  predicting  the  performance  of 
combustors  have  been  based  on  finite  difference  and  stirred  reactor  methods.  Du 
to  the  limitations  of  present  computers  finite  difference  methods  cannot  be 
extended  to  completely  include  fuel  spray  dynamics  and  realistic  chemical  kineti 
This  difficulty  is  overcome  by  using  the  computec  flow  patterns  to  define  a net- 
work  of  interconnected  stirred  and  plug-flow  reactors.  The  detailed  kinetic 
scheme  presently  consists  of.  13  species  undergoing  18  reactions  to  represent  the 
combustion  of  hydrocarbon  fuels  such  as  keroslne.  A model  of  fuel  evaporation 
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is  incorporated  which  assumes  the  fuel'  spray  to  be  composed  of  21  x 12  pun  size 
ranges,  evaporation  is  calculated  using  a forced  convection  model.  The  mixing  of 
fuel  vapour' and  air  is  modelled  using  a micro-mixing  parameter,  Tcp»  based  on  ' 
turbulence  dissipation  rates.  The  over-all  method  thus  combines  3-D  fluid 
dynamics,  turbulent  mixing,  evaporation  and  chemical  kinetics.  The  model  has  bee 
verified  by  experiments  which  show  that  flow  velocity  profiles,  chemical  species 
(CO, NO  etc)  can  be  successfully  computed. 
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